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EEFLECTION AND INTERFERENCE OF ELECTROMAGNETIC WAVES 

IN INHOMOGENEOUS MEDIA 

F. E. Geiger* and H. L. Kyle 
Laboratory for Meteorology and Earth Sciences 
NASA-Goddard Space Flight Center 
Greenbelt, Maryland 20771 


ABSTRACT 


Solutions were obtained of the wave equation for a plane horizontally polarized 
electro- magnetic wave incident on a semi infinite two dimensional inhomogeneous 
medium. The complex dielectric constant of the medium was assumed to be 
e (z) = e' (0) (1 + az) - i t" (0) (1 + bz), where a, and b are positive constants. 

Two problems were considered: an inhomogeneous half space, and an inho- 
mogeneous layer of arbitrary thickness for 0 £ z 5 Zi» contiguous to an inho- 
mogeneous half space for Zj < z < «. The dielectric constant e(z) is dicon- 
tinuous at z = 0 but continuous at z = z^, the gradient, grad e (z), is discontinuous 
at z = 0, Zj. Solutions of the wave equation were obtained in terms of Hankel 
functions with complex arguments. Numerical calculations were made of the 
reflection coefficient R at the interface of the homogeneous medium, z = 0, as 
a function of a, b, e‘(0), e" (0), by programming the Hankel functions and the 
expressions for R on a 360/91 computer. The startling results are that the 
reflection coefficient for a complex dielectric constant with gradient, e‘(0) a - 
i e"(0) b, can be less than that of the same medium with zero gradient, i.e, 
e(z) = e' (0) - i e"(0). The physical explanation of this behavior is given in 
terms of the interference of the coherent scattering reflections which take place 
thruout the inhomogeneous medium. 


*Prescnt address: 1301 Delaware Avenue, Washington, D, C. 20024. 
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REFLECTION AND INTERFERENCE OF ELECTROMAGNETIC WAVES 

IN INHOMOGENEOUS MEDIA 


I. INTRODUCTION 

The problem of a plane electromagnetic wave incident on a semi infinite inho- 
mogeneous medium has a long and interesting history. The problem is of two 
kinds, the purely analytical one of solving the wave equation with a variable 
coefficient describing the inhomogeneity of the medium and the physical inter- 
pretation of the results. 

Exact solutions have been obtained In a small number of cases, and various 
approximations which place a limit on the variation of the dielectric constant 
over the wavelength. The physical interpretation of the results would appear to 
be fairly clear, early doubts as to the presence of reflections in the inhomogeneous 
medium were dispelled by Wallot but primarily by Bremmer^, Similarly "re- 
flections" at discontinuities of gradients, and the question of unique resolution 
of the wave field into incident and reflected waves have been clarified^. But the 
fidl implications of the effect of internal reflections in the inhomogeneous 
medium seem not to have been appreciated. 

Microwave measurements of reflectivities of soils, and soils with varying 
amounts of moisture as a function of depth, from satellite and aircraft platforms 
made it necessary to get both insight into the physical nature of the problem, 
and actual figures of reflection coefficients. We therefore investigated the 
solution of the wave equation for semi infinite and layered media with linearly 
increasing dielectric constant in some detail. The solution of the problem for 
a real dielectric constant was obtained by Gans"* , and Wallot 

Wallot obtained a closed solution in terms of Hankel functions for an inhomoge- 
neous transition layer between two homogeneous media. The real dielectric 
constant is assumed to increase continuously from the free space value to 
the value e ^ of the second homogeneous medium. As the transition layer is 
reduced to zero thickness the Fresnel reflection coefficient is found to increase 
monotonically from zero to the limiting value corresponding to the Fresnel co- 
efficient between two homogeneous media of dielectric constants £q and e 
respectively. If we suppose that the transition layer is replaced by a single 
semi infinite inhomogeneous layer, Wallot's results clearly indicate that the 
reflection coefficient at the interface increases monotonically to the limiting 
value of 1 as the gradient approaches infinity. 

In our formulation of the problem we shall use the mathematical formalism of 
Brekhovsklkh® as the most suitable for our purposes. 
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Brekhovskikh calculates the reflection coefficient for a layered, and semi 
infinite inhomogeneous medium. Like Wallot, he uses a real dielectric constant 
and continuity of e (z) at the interface between homogeneous and inhomogeneous 
layer. In our calculation, we will change the physical characteristics of the 
inhomogeneous layer in three important respects, a) the dielectric constant is 
discontinuous at the first interface between homogeneous and inhomogeneous 
layer or inhomogeneous semi infinite medium, b) the dielectric constant is 
complex, c) the real and imaginary parts of the dielectric constant vary linearly 
with z, but independently of each other. As shown above, Wallots results indi- 
cate then that the reflection coefficient for an inhomogeneous half space with a 
complex dielectric constant will only differ from those with a real constant in 
detail but not in essentials. In other words we expect the absolute reflection 
coefficient to increase monotonically to the limiting value 1 as grad e (z) in- 
creases without limit. 

In fact the reflection coefficient behaves in a quite imexpected fashion depending 
on the individual gradients of e' (z) and e"(z). The reflection coefficient is no 
longer a monotonically increasing function of e'(z) and e" (z)*. We will show 
that the behavior of the reflection coefficient can be explained by interference 
phenomena in the inhomogeneous medium, and that Wallot ’s conclusions apply 
only to a special case. 

Early investigators were handicapped in the interpretation of their results by 
the unwieldiness of their results for the reflection coefficient. Hankel functions 
had only been tabulated for integral and half integral orders. Consequently only 
a small number of calcxilations were made under special limiting conditions, 
i.e, very large or very small arguments of the Hankel functions. We have made 
machine calculations of Hankel functions for complex arguments of order 1/3 
and 2 / 3 , and programmed our reflection coefficient calculations on an IBM 360/91 
computer. 


II. THE WAVE EQUATION AND WAVE FIELDS 
FOR AN INHOMOGENEOUS MEDIUM 

Consider a plane wave incident from free space on a semi infinite inhomoge- 
neous medium whose dielectric properties are a function of z only. (See Fig- 
ure 1.) The wave is polarized in the y-direction, the electric fields in the 


We were not able to prove analytically the monotonic increase of the absolute value of the 
Fresnel reflection coefficient with increasing e* and e" of the dielectric constant e' ^ i 6", but 
computer runs of |r | as a function of e’ and 6*' tend to show this. 
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X- and z-directions are zero, i.e, the incident field is horizontally polarized.^ 
The inhomogeneous halfspace may be separated by an inhomogeneous transition 
layer between 0 £ z 5 z j from free space. Both inhomogeneous media have a 
dielectric constant which is a linear complex function of z. Designating free 
space, the inhomogeneous layer, and the inhomogeneous halfspace as media I, 

0 II, and III, respectively, e (z) for medium II will be,* 


e(z) = e'(z) - ie"(z) = e'(0) (1 + az) - i e"(0) ( 1 + bz) 

0 £ z £ Zj, a £ 0, b £ 0, 


and in medium in, 

€(z) = e'(0) {(1 + azj) + c(z - Zj)) -ie"(0) {( 1 + bz^) + d (z - z i)) , 

It will bG seen that 6 (z) is a continuous function, but the gradient, Ve (z), is 
discontinuous at z = z^. 

The electric field satisfies the general wave equation® 

+ Ey(k^(z) - sin^(9o) = 0 (^) 

where , 

k^(z) ~ o)^/^(j£(z) 0 £ z < oo, media II and III, (2) 

and, ^0 > atigl© of incidence of plane wave. 


*MKS units and a positive time factor C}qj(i<wt) are used thruout. See References 6 and 7. 



In medium I, k^(z) reduces to the free space wave number e^, where, 


oi - 2 tt V , 


the circular frequency, sec" ^ , 


IJ,Q = 1.257 X 10-^ Hm'i, 


the permeability of free space. 


= 8.854 X 10-^2 Fm-i, 

the dielectric constant of free space. 

The solution of Eq. (1) in medium I is given by the incident and reflected wave- 
field, 


“ exp (- ik(j cos z) + R exp (ik(j cos ^0 z), (3) 

where, R is the reflection coefficient at z = 0. In medium II, we have the wave 
equation. 


+ kg2 [e’(0)/(l + az)/e^ - ie"(0)(l + bz)/eQ - sin^ ^ 0 


where, 


(4) 


^0^ - "Voeo- 


the free space wave number 
and, 
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k^(z) 


(e' (0) (1 + a z) - i e"(0) (1 + bz)) , 


( 5 ) 


e(0) = e'(0) - i€"(0), 

the dielectric constant of the inhomogeneous medium at the interface z = 0, 

We introduce the variable, 

^ = €‘(0) (1 + az)/6(, - ie"(0) (1 + bz)/eo - sin^^o, (6) 

and rewrite Eq, (4), 

B2Eya^2 + ^k„/(ae-(0)/£o * i b e"(0) /^o)) ' ^ = 0. (?) 

We make a second transformation. 


w 


2/3 


ko \ 

ae'(0)/eo - ihe"(0)/€^)^ 


( 8 ) 


writing the wave equation in terms of the new variable w, 

+ 1/3 w-‘ 3Ey/3w + Ey = 0. (9) 

The solution of Eq, (9) can be written most conveniently for our purpose in 
terms of Hankel functions® of the first and second kind of order 1/3, 


Ey = (w) + B (w), ^ 0 £ z £ Zj, (10) 

Similarly we have for medium III, 
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( 11 ) 


[e'(0) ((1 + azj) + c (z - Zj)) /e^. 


- i e"(0) ((1 + hz^) + d (z - zi)) /e^ - sin^ = 


Making the transformations as above. 


V - e'(0) ((1 + azj) + c (z - z^)) /e„ 


( 12 ) 


- i e"(0) ((1 + bzi) + d (z - z^)) /e^ - sin^ 5(j, 


V = 2 / 3 \ 


c e'(0)/eo - ide'(0)/€. 


, 3/2 


(13) 


= Cv'/"H,/n)(v) + Dv1/3Hi/3<2)(v), (14) 

where C and D are constants as are A and B in Eq. (10). 

Ha. Determination of the Constants A, B, C, D and the Reflection Coefficient R 

The electric field must remain finite as z goes to infinity. From Eqs. (12) 
and (13) we have, 

z"- I = 2/3 ko z^/2 i ^ 

or V - as z -> «>, and v will be in the fourth quadrant, 

(- 1/2 7T < arg V < 0) for c > 0, d > 0, > 1, 
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hence «•, and -* 0 as v -* «>, and we must have^° C = 0. And 

we have, 


Ey ^ (v), Zj < z < «. (15) 


The continuity conditions for and SE^/Bza at z = give the following 
equations, 


A (w) + B (w) - D 2) 2 - z^, 

(16) 



vi 


- V 


Z-Zj 


W, 


^ (Aw^/3H_2/3(D(w) + Bwl/3H_2/i2)(w)) - ^ Dv1/3h_2/3(2)(v), 

(17) 


where we have used the relation,^* 




W^'^3h_2/( 1).(2) (w). 


(18) 


From Eqs. (6), (8), (12), and (13), 


^ li . k .1/2 ^ ±!Z . k 7,1/2 

dz ’ dn dz *^0^ 


( 19 ) 


and 


dw _ dv 

d z d z ’ ^ ^ 
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Solving Eqs. (16), and (17) for D, 

D = (Aw1/3H,/3(1)(w^) + Bw,i/3Hj^3(2)(wi))/(v^3H,/P)(v,)). (20) 

we find for A/B from Eq. (16), 


A/B 


H-2/3 ('''l 


( 21 ) 


The continuity conditions at the interface z = 0 give from Eqs, (3) and (10), 


( 1 + R) = A ^Hw) + B Hj/i (w) 


( 22 ) 


z-0 


and 


w„ - w 


Z “0 


and 


iko(R - 1) cos ^ 


z- 0 
w=w« 


dw 

dz 


z- 0 


- ko p(0)/€o - i€"(0)/€, - 


1/2 


(23) 


(24) 


From Eqs. (23) and (24) we easily find for the reflection coefficient at z =0, 


R = 


- i cos 9 q - ^e'(0)/eQ - ie"(0)/€Q 


sin 


• G 


0 / 


i cos 0(J + ^£'(0 )/€j) - ie"(0)/£o - * 


(25) 
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where, 


(A/B) ^ 

(A/B) 


Equation (25) looks very much like the Fresnel equation for reflection of a plane 
wave from a semi infinite medium with dielectric constant €'{0) - i e"(0), with a 
factor G (Wg , Wj, Vj) in the numerator and denominator to account for the in- 
homogeneity. In fact we can show that Eq. (25) reduces to the Fresnel equation 
if the arguments of the Hankel functions are allowed to go to infinity, i.e, a — b 
= c = d - 0, then Iw^l - «, | wj - and jvj - « . In other words the in- 
homogeneity has been removed. 

Using the asymptotic expansion of Hankel functions^ ^ for large arguments, we 
find, 


= ( 2 /t 7 exp (i z - i7r5/12) (l - i • 5/(72z)), (27) 



(2/t7 z) 


(2/77 Z) 


(2/t7 z) 


(2/77 z) 


where all arguments of z are assumed to be in the 4th quadrant. After a great 
deal of algebraic manipulation we have. 


A/B = (1/12) exp(-i2wj + i77/3) ^1/w^ - 1/vj + 0 (w v (31) 

and 

- i (l - i/(6W(,)} + (1/12) exp (2i(w„ - Wj)) (1/w^ - 1/vj) 

G - ^ ^ (32) 

1 - (i/12) exp (2i(wQ - wj) (l/w^ - 1/vj) 
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In the limit as 1 Wq I , I Wj I , | | - ~ , 

lim G “ - i , 



ao 


and Eq. (25) reduces to the Fresnel equation, 


cos ^0 ■ (c'(0)/eo - ie"( 0 )/eo - sin^-9J 

R = . (33) 

cos 6>o + ^€'(0)/€(, - ie"(0)/€o - 


lib. Limiting Expressions for the Reflection Coefficient 


The double layer expression for the reflection coefficient R, Eq, (25), can be 
reduced to a single inhomogeneous layer (i.e, inhomogeneous half space) ex- 
pression by settii^ c = a, and d = b, then remains unchanged, and Wj = v^. 
Equation (21) reduces to zero and G becomes. 


G = H_2/i2)(wo)/Hi/3<2)(wo). • (34) 

If the first layer, (0 5 z < z^), is inhomogeneous, and the half space, (Zj < z < «.), 
homogeneous, c = d = 0,Vj and we have the following expressions for A/b 
and G, 


i ^l/i ("'l) ^-2/3 

H-2/3^nWi) + iHi/3(2)(Wj) 


(35) 
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No useful limiting expression is found if the first layer is homogeneous and the 
second an inhomogeneous half space. In this case a = b = 0, andw^ approai 
infinity. This problem has to be solved by considering the specific boimdary 
value problem from the beginning, and solving the equations for the reflection 
coefficient, 

lie. Machine Calculation of the Reflection Coefficient 

The expression for the reflection for an inhomogeneous medium is so complex 
that no conclusions can be drawn about the general behavior of the coefficient. 
Equations (25), (26), and (21) were therefore programmed for an IBM 360 com- 
puter. The program included also direct calc\ilation of the Hankel fxmetions of 
the first and second kind of orders 1/3 and 2/3. Tabulations of these functions 
are not available except for modified Hankel functions with complex arguments 
and absolute values | z] 56 . This is much too restrictive a range for the prob 
lem under consideration, and for a general exploration of the behavior of R. 

We therefore wrote our own straightforward computer program using the fol- 
lowing equations to calculate the Hankel functions for an arbitrary complex 
argument.^^ For |z| < 7.5 we used, 


= 

i CSC (V77) [exp (- ITTV) J^(z) - J-„(z)] 

(37 

H/2)(z) = 

- i CSC (V77-) [exp (iTTV) (z) - (z)] , 

(38 


where the Bessel functions were calculated using the ascending series. 


Jv(^) 


(1] " y (- 1/4 zT 

\2/ ^ m ! r (v- + m + 1) 

m^O 


(39 


The series was summed in double precision arithmetic to yield J (z) to 6 1/2 
significant figures. For 1 z 1 > 7.5 we used terms throi^h y in the asymptotic 
expansion of the Hankel functions: 


- '^ 1 / {jr z) [P (v, z) + iQ(v, z)j exp(iX) (40 
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" y'2/(77- z) [p (v, z) - iQ (i^, z)] exp(-ix), 


(41) 


where, 


X - z - (l/2v + 1/4) 77-, /^ - 


I 


P(v. z) - 


k = 0 


(u, 2 k) 

(2z)2k 


Q (i', z) - 


E (- 1)" 

k = 0 


(v, 2k + 1) 
(2z)2k + l 


and, 


(v,k) 


12) (^- 32) .... (m- (2k- 1)2) 
22k k ! 


(i^.O) - 1. 


This prescription yields Hank el functions accurate to six significant figures 
for I z I < 7.5 or | z | > 10, and five significant figures for 7.5 < 1 z | < 10. 
Greater accuracy can be obtained by retaining more terms in the expansions. 

III. RESULTS OF COMPUTER CALCULATION 
OF THE REFLECTION COEFFICIENT 

Machine calculations were made of the reflection coefficient for single and 
double layer inhomogeneous media at several microwave frequencies. The 
calculations on the single layer semi infinite inhomogeneous medium lend them- 
selves more readily to interpretation and understanding of the underlying 
physical processes, and will be discussed first. 

Figures 2,3, and 4 show a plot of the reflection coefficient for an inhomoge- 
neous half space and horizontal polarization. The dielectric constant at the 
vacuum- dielectric interface was chosen to correspond to soil with roughly 
2% moisture, and alternately with 15% moisture by weight. Moisture content 
was assumed to increase as a function of depth (distance from interface) re- 
sulting in a hypothetical linear increase of e'(z) and e’'(z), or constant gradient 
d(z)/dz = e' (0) a - i e''(0)b. The absolute amplitude reflection coefficient was 
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plotted as a function of the constant b, which controls the gradient of the 
imaginary part of the dielectric constant, for various values of the constant a. 

Intuitively one expects | R | to increase from the value determined by the Fresnel 
equation for e = e'(0) - i e"(0) with increasing a and/or b. Although there is an 
immediate increase in |R I for a = b > 0, further increase in b results in a 
steady decrease of 1 R | until a minimum is reached. We find the surprising 
fact that for a broad range of values of a and b the reflection coefficient at the 
interface drops considerably below the value of the Fresnel coefficient of a 
semi infinite homogeneous medium with dielectric constant (e'{0) - ie"(0)) , i.e. 
the dielectric constant of the inhomogeneous at the vacuum- dielectric interface. 
On the other hand if b = 0, but not necessarily €"{z) the reflection coefficient 
increases steadily with increasing a from the aforementioned value of the re- 
flection coefficient. This behaviour is shown in Figure 5, and is entirely in 
accord with one's physical intuition. 

One is forced to conclude, since the absolute value of the Fresnel coefficient is 
a monotonically increasir^ function of e'(z) and e”(z)*, that the drop in the re- 
flection coefficient of the inhomogeneous medium is caused by interference, 
Bremmer ^ explored the physical interpretation of the wave field in inhomoge- 
neous media, and was able to demonstrate the presence of reflected waves at 
all points of the inhomogeneous medium. We will try to demonstrate in a sub- 
sequent section that the phase relationship of the reflections in the inhomoge- 
neous medium is such as to cause interferences and a reduction in the amplitude 
of the wave emerging from the inhomogeneous halfspace. 

A limited number of calculations were also made for a two-layered inhomogenous 
medium at three microwave frequencies, 37,0, 19.35, and 5.0 Gc/ sec. The first 
inhomogeneous layer is of finite thickness and is followed by a semi infinite in- 
homogeneous medium. The dielectric constant is continuous at the interface 
between the two inhomogeneous media, but there is a discontinuity in the gradient 
of the dielectric constant, i.e, 3e/9z = e’(0)a - ie"(0)b for 0 < z 5 Zj, and 
3e/3z = e'(0)c - ie"(0)d for Zj < z < »o, where z^ is the distance from the 
free space-dielectric interface to the interface between the inhomogeneous 
media. 

Figure 6 shows the results obtained for a two layer model for two values of z 
However the dielectric constant gradients and the dielectric constant e (0) = 
e'(0) - i e"(0) at z = 0 were the same in both cases. For comparison we calcu- 
lated the reflection coefficient as a fxmction of frequency for a simple semi 
infinite inhomogeneous medium whose dielectric constant e (z) is the same as 


*See footnote on p. 3- 
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that of the first layer of the above composite (two layer) model. There is a 
striking difference in the behavior of the reflection coefficient for the three 
cases. The single layer reflection coefficient is relatively insensitive to changes 
in frequency, the double layer, on the other hand, shows relatively large fluctua- 
tions in the coefficient when is changed from 0.5 x to 0,25 x 10-^ meter*. 
We conclude, therefore, that the discontinuous change in the gradient Be/B z at 
z = Zj gives rise to strong reflections, which emerge at z = 0 and interferes 
with the reflection at the air-dielectric interface, A change in z ^ of 0.25 cm 
corresponds very roughly to a quarter wave length for 19 Gc/ s and an average 
dielectric constant of e = 3,0 x 10"^^. Thus the emergent reflection will have 
suffered an approximate phase shift of tt, and will be in phase with the reflection 
at the interface. This reasonii^ assumes that a definite reflected wave can be 
assumed at the interface z = z However the wavefield in the first inhomoge- 
neous layer cannot be separated into incident and reflected waves so that a re- 
flection coefficient at z = Zj has no meaning and cannot be defined. 


IV. INTERFERENCE EFFECTS IN INHOMOGENEOUS MEDIA 

The approximate solution of the wave equation (see Eq. (1) ) for an inhomogeneous 
medium extending from -« to +» with a slowly varying wave number k (z) is the 
well known WKB solution*^ , 


- i k (z) d 
'o 

I 

where C^ and Cj are constants. There are no reflections to this order of 
approximation, and the two waves in the positive and negative directions, re- 
spectively, are independent of each other. Bremmer^® has shown that Eq. (42) 
represents only the zero order term in a series solution, and higher order 
approximations can be interpreted as reflected waves. Thus the zero order 
wave 





propagating in the positive z-direction from -«>, produces a reflection at each 


^Suggested by T. Wilheit. 
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point in the medium given by 


cp> 




f (z) exp 




Conversely, the zero order wave 


Cf > 



produces a first order reflected wave propagating in the positive z-direction, 
etc. It can easily be shown that the zero order waves propagate without reflec- 
tion. The power flow calculated from the Poyting vector 


S = (EyXHJ 

for the "upgoing" wave 



is a constant, i.e, undiminished by reflections. In order to demonstrate the 
presence of interference effects in a semi infinite medium, as apparently 
manifested by the drop in the reflection coefficient at the vacuum-dielectric 
interface, one would have to show simply the presence of emerging scattering 
reflections such as C } . etc. These obviously coherent reflections 

are then capable of interference. Their exact phase relationship will be deter- 
mined solely by the dielectric constant e (z) of the semi infinite inhomogeneous 
medium. 

However, for a boundary value problem such as ours it is no longer possible to 
demonstrate in a perfectly general way the reflected waves in the inhomogeneou 
half space. The WKB solution of the wave equation for a plane wave incident on 
a semi infinite inhomogeneous space is. 
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z > 0. 


E 


y 


Cl (k(z))~^^^ exp 



Using Bremmer's approach in the preceding paragraph for an infinite medium 
(-<» < z < +«) and applying it to the semi infinite case (0 < z < ») turns out to 
be not very fruitful. It seems impossible to arrive at higher order approxima- 
tions to Ey, which could be identified as secondary, . . . etc. waves originating 
from reflection losses of a primary, . . . etc. wave. A more direct approach 
suggested by Bremmer's exact solution of a plane wave reflection from a semi 
infinite medium^ ^ with k(z) a l/z consists in examining the reflection coefficient 
(see Eqs. (25) and (34) ) directly. 

Using the as 5 rmptotic expansion of 2) and H w 1 » 1 of Eqs. (40), 

and (41) one finds for G, (see Eq. (34) ), 


^ ~ ^-2/3^^ (Wo) - exp(- iTT/ 2 ) . 

(1 + 1 bj/w(, - b2/W(2 . . . ) 

Using the binomial expansion, we have, 

G exp (-177/2) [l - (aj + bi)/wo + (aj + bj ~ ajbj - + 0 (wg^)) 

where, 

= 7/72 

bj ^ 5/72 

32 = 455/10368 
bj ^ 385/10368. 

From Eq. (25) we have for the reflection coefficient, for normal incidence. 
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1 “ (^r (0))^^^ “ i (3i + t>i)/wo + (32 + bj - a^bj - hi)/-WQ . . .) 

1 + (0)) (l - i (aj bi)/wo + etc. . . .) 

where 

K(0))'"' - (e’(0)/€o - ie"(0)/6o)'''. 

We can rewrite R by again using the binomial expansion and assuming 
(er(O))^''^/ (l + ^ t>i)/Wo| , 

and 

|(e,(0))^''V (l + e,(0))^^^(a2 + bj - a^bj - \ « 1- 

Retaining only second order terms in Wj^, and after much algebraic manipulation, 
we find, 

R = (l - (l + (i/6) 2(€.(0))'"V(1-.,(0))/»„ 

- <l/»o“) • ((0.083)2. (€,(0))‘"V(l - £, (0)) HO. 028)2 .e,(0)/ (l - £“(0))) . . . j . (43) 

The first term in the brackets clearly corresponds to the WKB approximation. 

It is simply the reflection coefficient for a homogeneous semi infinite medium 
with relative dielectric constant 

e,(0) - €-(0)/6o - i e"(0)/6^, 

and as shown above there are no contributions to the scattering reflections from 
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the inhomogeneous medium in the WKB approximation. One can also prove the 
absence of scattering reflections by solving the boundary value problem of the 
semi infinite inhomogeneous medium in the WKB approximation*. The result, 
as anticipated is simply the reflection coefficient for a homogeneous medium of 
dielectric constant e^. (0), The second and third terms of Eq. (43) are complex 
constants whose magnitude and argument represent the amplitude and phase, 
respectively, of the scattering reflections emerging at the interface and propa- 
gating in the negative z-direction. In other words the first secondary reflection 
due to the primary progressive wave, (see Eq. (42)), 


E 

y 


Cl (k(z)) exp 



(44) 


in the inhomogeneous medium is given by. 


l+(e^(0))‘^" (l-e^(O)) 



exp (i k^z). 


(45) 


These secondary reflections combine destructively or constructively depending 
on the phase relative to the primary reflection at the interface due to the WKB 
term. See Eq, (43). 

We have thus demonstrated the presence of secondary plane wave reflections in 
the homogeneous subspace -<» < z < 0 and by inference secondary reflections in 
the inhomogeneous subspace 0 < z < «>. 


V. THE RADIATION FIELD IN THE SEMI INFINITE 
INHOMOGENEOUS MEDIUM 

In the WKB approximation (see Eq. (42) ) and the Bremmer solution for an 
infinite inhomogeneous medium waves in the positive and negative direction can 


*See Appendix. 

t It may be shown that the exact solution for (see Eq. (10)) in terms of H i/ 3 ^^^(w) reduces to 
a first approximation to the WKB solution, Eq. (44). See Reference 18. 
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still be identified, although some may be independent of each other as shown in 
the preceding section. The exact solution of the wave equation for the specific 
case of a linearly increasing dielectric constant in a semi infinite space as 
given by Eq. (10), 


Ey = Hj/<2)(w), 0 < z < cc, 


no longer permits a clear cut distinction between reflected and incident waves, 
or even waves travelling in positive and negative directions independently of 
each other. 

In order to get a more physical picture of the wave field in the semi infinite 
space, and in addition the powerflow into the medium the complex Poynting 
vector S* was machine calculated from 

iReS;| - lEe|(E>:H)l = <^6) 


From Eq, (10) and from Maxwell's equation for B E^/Bz, we find 

for normal incidence, = 0, 


Then, 


Re S’ 


1 

Re 2 




where the sign denotes the complex conjugate. 

If the medium is not lossy machine calculations show that Re S* is a constant 
for all z >0^, The discussion of internal reflections throughout the inhomoge- 
neous medium of the preceding section would have led one to guess at a constantly 


^ On the basis of a number of computer runs for very large a’s (see Eq. (5)) (the gradient of the 
dielectric constant is proportional to a), and distances up to 5.0 cm from the interface. 


19 



decreasing Re S* as z - <«, But some reflection shows that in fact a constant 
Poynting vector must be correct. If Re S* - 0 as z - °o, and no energy is dis- 
sipated in the medium all the reflected energy in the inhomogeneous medium 
would eventually have to emerge at the interface with the homogeneous medium, 
resulting in a reflection coefficient of 1. This is absurd. And contrary to in- 
tuition the energy flux in the positive z-direction stays constant as in a homoge- 
neous lossless dielectric medium,^ 

Results of calculations of Re S* when the medium is lossy, e" (z) / 0 , are 
shown in Figure 7 for three different dielectric constant gradients. Curves I 
and II of Figure 7 correspond roughly to an Arizona loamy soil with approxi- 
mately 15% and 2% moisture by weight at the interface. It is clear from the 
figure that the functional dependance of the powerflow on penetration distance 
is not the same for all curves. Curve I happens to be very closely exponential. 
Curve III roughly follows an exp (-ax^) decay. 

Losses in the inhomogeneous medium may be calculated directly from the 
expression^ 


-If*'" 

Loss - J Ey EyCr dz, 


(49) 


where cr(z) is the conductivity of the medium, which may be expressed in terms 
of the dielectric constant cr(z) = f "(z) Using Eq. (10) for E^, and the ex- 
pression for cr(z). 


Loss = I 6;e"(0) JJ(1 + az) !b|^ ■ iHj/j 2) (w)l^ dz . (50) 


This integral was machine calculated by Simpson's rule. The calculations show 
the normalized L to be exactly equal to (1 - Re S* ). Additional proof that any 
decrease in the Poynting flux Re S * in the inhomogeneous medium is solely due 
to the conduction (heating) losses and not any scatterii^ reflection losses. 


^This result was anticipated by Waliot, who gave an analytic proof of the constancy of the energy 
flux in the transition region between two homogeneous media for a real dielectric constant. See 
reference 20. 
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VI. CONCLUSIONS 


We have shown that the reflection coefficient of a semi infinite inhomogeneous 
medium with constant positive gradient de (z)/dz = e’ (0) a - i e" (0) b is not 
necessarily a monotonically increasing fxmction for increasing a and b. Machine 
calculations of the reflection coefficient I R I , show the coefficient to decrease 
with increasing b until a minimum is reached. On the other hand | R | seems to 
increase monotonically with a for a given b as far as can be judged from the 
available machine calculations. 

We demonstrated the presence of scatterii^ reflections in the free space 
-«. < z < 0, emerging from the inhomogeneous halfspace, in analogy with 
Bremmer's results on infinitely and semi infinitely extensive media. But in 
addition we drew the important conclusion that these scattering reflections, 
which are coherent, interfere to produce the calculated drop in the reflection 
coefficient despite increasing gradients of the dielectric constant. As a conse- 
quence we find the startling result that a semi infinite inhomogeneous medium 
with a steep dielectric constant gradient may have a much smaller reflection 
coefficient than the same dielectric half space with zero gradient, i.e. a = b = 0, 
and e(z) = e'(0) - ie" (0). 
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APPENDIX 


We shall consider the reflection of a plane horizontally polarized wave, at 
normal incidence, from an inhomogeneous semi infinite medium in the WKB 
approximation. The dielectric properties of the inhomogeneous medium are 
characterized by the propagation constant (z) = ^ (z). (See Eq, (21),) 

Then from Eqs, (42) and (44), we have for the progressive wave in the inho- 
mogeneous medium. 


E 


y 


C (k (z)) exp 



(Al) 


and in the homogeneous medium as before from Eq. (3), 


Ey - exp (- ik(, z) + R e:iq3 (+ ik(, z). (A2) 

From the continuity conditions at z = 0, we find from Eqs. (Al) and (A2), 


1 + R = C/(k(0))^^^ 


(A3) 


^ C dk(z)/dz _ C k(0) 

iko(k(0))'/'* (k(0))'/'k„ 


(A4) 


In the WKB approximation (dk (z)/dz)/k(z)^''2 ~ O, then solving Eqs, (A3), 
and (A4) for R, we find. 


ko - (k(0)) 

R " 

ko + (k(0)) 


(A5) 


Equation (A5) is the reflection coefficient for a homogeneous medium with 
dielectric constant e (0), Consequently, in the WKB approximation there are 
no scattering reflections in the inhomogeneous medium. 
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Figure 1. Geometry of a plane wave incident from free space on the x-y face of a 
semi infinite inhomogeneous medium. The plane of incidence is the x-z plane and 
the wave is polarized in the y*dlrectlon. There is an inhomogeneous transition 
layer between 0 S 2 1 z^and both inhomogeneous media have a dielectric con- 
stant which is a linear complex function of z. Free space, the inhomogeneous 
layer and the inhomogeneous half space are designated as media l,JI/ and III 
respectively. The angle of incidence, Sq, between the wave vector and the 
normal to the plane is taken to be zero in all cases: The permiability and dielec- 
tric constants of free space are respectively and £q. 
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ABSOLUTE AMPLITUDE REFLECTION COEFF (PERCENT) 



Figure 2. Absolute amplitude reflection coefficient for a semi infinite inhomogeneous medium 
with dielectric constant linearly dependent on z. No boundary layer; horizontal polarization; 
angle of incidence, 0^; frequency, 19.35 Gc/sec. 
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ABSOLUTE AMPLITUDE REFLECTION COEFF (PERCENT) 



Figure 3. Absolute amplitude reflection coefficient for a semi infinite inhomogeneous medium 
with dielectric constant linearly dependent on z. No boundary layer; horizontal polarization; 
angle of incidence, O”; frequency, 19.35 Gc/sec. 
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ABSOLUTE AMPLITUDE REFLECTION COEFF (PERCENT) 


V = 5.0 Gc/sec 


19.35 Gc/sec 


HORIZONTAL POLARIZATION 
ANGLE OF INCIDENCE 0® 

DIELECTRIC CONSTANT €{Z) = €' (o)(1 + aZ) • i e" (oHI + bZ) 
< ' (o) = 3.0 X 10 FM ■’ , £ " (o) = 0.2 x 10 ” FM ’ 


10,000 20.000 30,000 40,000 

b. METER •’ 


50,000 


60,000 7( 


Figure 4. Absolute amplitude reflection coefficient for a semi infinite inhomogeneous medium 
with dielectric constant linearly dependent on z. No boundary loyer; horizontal polarization; 
angle of incidence, 0 ; frequency, 19.35, and 5.0 Gc/sec. 




ABSOLUTE AMPLITUDE REFLECTION COEFF iPERCENT) 



a, METER-’ 


Figure 5. Absolute amplitude reflection coefficient for a semi infinite inhomogeneous medium, 
with dielectric constant linearly dependent on z, as a function of the parameter a. No boundary 
layer; horizontal polarization, angle of incidence, 0°; frequency, 19.35 Gc/sec. 
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UUt HtLtUI lUh 


Z, =0.25 CM 


SEMI INFINITE MEDIUM 


Z, =0.50 CM 


TWO LAYER INHOMOGENEOUS MEDIUM 
HORIZONTAL POLARIZATION, = 0° 

«(Z)= e'(o) [(1 +aZ, )+c(Z-Z,)).i€"(o) 1(1 +bZ,) +d{Z-Z,)J 

€* |o) = 2.8x 10 " FM \ €" (o) = 0.078 X 10 ’’ FM ’ 
a = 80, b = 1800, c = 434, d = 15180 

SEMI INFINITE MEDIUM e(Z)=€'|o)(1 +aZ)-ie"(o) (1 +bZ) 


FREQUENCY, v , Qc/sec 

Figure 6. Absolute amplituc'e reflection coefficient for a two layer inhomogeneous medium os o 

function of frequency. Dielectric constant of inhomogeneous media lineoHy dependent on z. 
Morizontaj polarization; angle of incidence, 0^. 






POWER NORMALIZED AT INTERFACE (FRACTION) 


LINEAR VARIATION OF DIELECTRIC CONSTANT WITH DEPTH. t(Z) = £'(Z) -i£"(Z) 


I £' (o) = 5.0 X lo ’"' , €" (o) * 1.0 X 10‘" 
a = 200, b = 200 

n €’ (o) = 3.0 X lO '"' ' , £" (o) = 0.2 X lO '’"' 
a = 50, b = 400 

(ADELANTA LOAM, ROUGH APPROXI 
NATION TO PROFILE H) 

HI £' (o) =3.0x 10 ’\ €"(o) = a2x10r’'’ 

a = 0, b « 7000 


Z, DEPTH OF PENETRATION, CM 


Figure 7. Power flow in semi infinite inhomogeneous medium with no boundary layer as a func 
fion of the depth of penetration. Frequency, 19.35 Gc/sec, normol incidence. 





